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From confining fields back to power corrections 

V.I. Zakharov a 

a Max-Planck Institut fur Physik, 
Fohringer Ring 6, 80805 Munchen, Germany. 

We overview the issue of the power corrections to the parton model, mostly within the context of QCD sum 
rules. There are a few sources of the power corrections responsible for various qualitative effects. The basic idea 
that the power corrections are related to confining fields seems to be realized in case of 1/Q 2 corrections which 
are dual to high orders of perturbation theory. The evidence comes from lattice studies of the confinement which 
allowed to identify the confining field configurations as monopoles and vortices. We comment also on the role of 
the condensate of dimension two, < (A^) 2 > min . 
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1. INTRODUCTION 

The issue of the power corrections within the 
context of the QCD sum rules P has been sub- 
ject of many reviews, see, e.g., [21314) . Many 
things have become common knowledge and there 
is no need to review them further. Moreover, we 
actually assume that the reader is familiar with 
sum rules, at least, in general terms. One of the 
main questions, however, has not been yet fully 
resolved. Namely, the power corrections were in- 
troduced first to parametrized the effect of 
confining fields. Nowadays, however, this issue 
is overshadowed by technical details of applica- 
tion of the sum rules. It might be still useful to 
evaluate the status of this issue. 

Actually, it is at this point that the situation 
now is radically different from the early years of 
the sum rules. Namely, there was little hope that 
the confining fields could be identified explicitly 
since confinement is a property of strong-coupling 
theory. Thus, the sum rules were intended to at 
least systematically parameterize effects of con- 
finement at short distances. Nowadays, lattice 
studies seem to reach the point when the knowl- 
edge of confining field configurations is quite pre- 
cise. The confining field configurations are the 
lattice monopoles and vortices, for a recent re- 
view see, e.g., [Sj. Thus, one can try to project 
these fields back to short distances and clarify 
what kind of power corrections they generate. We 
will argue that it is the so called quadratic cor- 



rections j(il7l,'{j which seem to be directly related 
to the confining fields [Sj. 

The issue of the quadratic corrections is some- 
what controversial. On the scientific side, there is 
a well defined reason for this. In terms of the Op- 
erator Product Expansion (OPE) the quadratic 
corrections are a part of the coefficient functions 
in front of the unit operator and are not related 
to a matrix element of a non-trivial operator. In 
other words, the quadratic corrections are actu- 
ally dual to high orders of perturbation theory. 
This is a unique feature of these confinement- 
related power corrections and the common cul- 
ture of working with the sum rules does not help 
much to appreciate this special feature of the 
quadratic corrections. 

On the other hand, this manifestation of dual- 
ity in pure Yang-Mills theories might help in un- 
derstanding dual formulations of Yang-Mills the- 
ories which is commonly believed to be a string 
theory, see, in particular, [§]. 

2. SUM RULES 

2.1. Classical channels 

In somewhat simplistic form, the sum rules 
read: 

ds exp(—s/M 2 )Rj(s) w (parton model) ■ 

•(1 + a 3 a s {M 2 ) + bj < G 2 > /M 4 ) (1) 

where M 2 is a parameter, M 2 ^> Aq CD , the 
coefficients a,j , bj are calculable perturbatively, 
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Rj is proportional to cross section induced by a 
hadronic current j (if, say, j is the electromag- 
netic current then Rj is proportional to the total 
cross section of e + e~-annihilation), a s (M 2 ) is the 
running coupling, < G 2 > is the gluon conden- 
sate. There are actually higher power corrections, 
say, proportional to M~ 6 which we omit here. 

It might worth emphasizing that the sum rules 
Q are far from being 'natural'. Indeed, we cut 
short the perturbative series by keeping only the 
hrst perturbative correction, aja s (M 2 ) and, on 
the other hand, keep the power corrections. How- 
ever, for very small a s any power of a s is much 
larger than any power correction since M~ 2 ~ 
eM~b /a s (M 2 )). 

The idea behind this breaking of the perturba- 
tive hierarchy was that the pure perturbative cor- 
rections are not signalling confinement. And the 
physics of confinement is encoded in the power 
corrections which become crucial at small, but 
not very small a s and are synchronized with de- 
viations from the parton model in the measurable 
cross section Rj. 

There is no a priori justification for such an 
assumption and one relies entirely on the phe- 
nomenological success to pursue the sum rules in 
the form QJ. In particular the sum rules in the 
p- and J/^-channels allowed to 'explain' the ap- 
proximate equality of the mass splittings: 

"V - m,//* ~ m p > - m p . (2) 

Although Eq (J2J) looks trivial it is actually diffi- 
cult to understand theoretically since relativisti- 
cally we should compare the masses squared, and 
these are very different in different channels. Sum 
rules based on the asymptotic freedom alone ^0] 
fail to explain (J2J completely. Explaining (J2J) was 
the first qualitative success of the sum rules QJ. 

2.2. Hierarchy of scales 

Sum rules Q have numerous successful appli- 
cations. There are, however, remarkable failures 
as well jlll4) . Namely, the most straightforward 
prediction of the sum rules is the value of M 2 rit at 
which violations of the parton model become siz- 
able. The numbers vary from channel to channel 
because the coefficients like bj (see ||TJ) depend 
on the channel. 



A careful (independent of the sum rules) phe- 
nomcnological analysis allowed to establish a 
kind of hierarchy of scales: 

« 4-(M c 2 rit ) p (3) 

{M cr it) giueball 0+ ~ ^iMcrit) ff 

The standard sum rules JD fail completely to re- 
produce this hierarchy. 

2.3. Direct instantons 

There is one quite obvious reason to modify the 
sum rules Indeed, none of the terms kept in 
would signal special role of the rj -channel. It 
is known since long |12| that one has to include in- 
stantons to explain shifting of the r\ -meson mass 
compared to other pseudoscalar mesons. 

The instanton contribution at large M 2 starts 
as a highly suppressed power correction ^2] to 
the right hand side of Eq : 

^instantons \ ^ j J • V / 

However, at values of M 2 which become sensitive 
to the resonance structure in the left-hand side of 
||TJ this correction should overcome other sources 
of the power corrections, at least in case of some 
channels, like the r] -channel. 

Detailed analysis of the instanton contribution 
relies actually on elaborated models of the vac- 
uum, like instanton-liquid model, for review see, 
e.g., 0]. There is overwhelming evidence in favor 
of the instantons in some channels. 

3. QUADRATIC CORRECTIONS 

3.1. Short-distance gluon mass 

Even upon inclusion of the instantons there re- 
main phenomenological problems with the sum 
rules. The example easiest to understand is 
provided by the heavy quark static potential at 
short distances. Strictly speaking, the potential 
cannot be treated within the (OPE). However, 
quite a straightforward extension of the OPE pre- 
dicts that the correction due to the soft non- 
perturbative fields is of order 

SV QQ (r) ~ A 3 QCD -r 2 (5) 

where r is the distance between the quarks and r 
is small, r ■ Aqcd <§C 1. 
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However, the lattice data are fitted by the so 
called Cornell potential, 



V QQ (r) 



const ■ a. 



(6) 



at all distances available. In particular, if one 
trusts this fit at short distances there is a linear in 
r correction to the Coulomb-like potential, 5V ~ 
a-r 0. 

The argumentation just presented is actually 
not so strong. The problem-quite usual one with 
the power corrections- is that one should con- 
fine oneself to the distances where the power cor- 
rection to the Coulomb-like potential is small. 
But then it is difficult to extract the hypothet- 
ical piece SV ~ a ■ r on the background of the 
perturbative corrections, see, in particular, |15| . 

The only way out, to my mind, is to analyze 
along the same lines a few channels or find quali- 
tative effects. In case of the quadratic corrections, 
unfortunately, there is no rigorous way to relate 
corrections in various channels. Indeed, so far we 
have not yet even identified the source of such 
corrections. 

An attempt to break the deadlock was made in 
ref. [H] by introducing a model to relate quadratic 
correction in various channels. The starting ob- 
servation is that to reproduce the potential © at 
short distances one can modify the gluon propa- 
gator at short distances by introducing the gluon 
mass: 



1/Q 2 -» 1/Q 2 + m 2 /Q 4 



(7) 



where Q is the Euclidean momentum carried by 
gluon and Q 2 » ^qcd- Ones replacement {7J| is 
postulated, the quadratic correction is calculable 
in other channels as well. 

It turns out that the model Q , using as an in- 
put the string tension a immediately reproduces 
both the overall scale and the hierarchy of the 
anomalous contributions J2J 6|. There are a few 
other non-trivial checks of the model. Note also 
that m 2 in Eq (JJJ is positive which implies that 
the short-distance gluon mass is in fact tachyonic. 
Later the prediction of the model were scrutinized 
further and so far the model works very well, 
for a review see the talk by S. Narison at this 
Conference ^7] 



3.2. Ultraviolet renormalon 

In view of the phenomenological success of the 
model J7J it is useful to clarify the theoretical 
status of the model. Historically, the power cor- 
rections were introduced first within the context 
of divergences of perturbative series, for review 
see, e.g., [T%] . 

Consider, for example, polarization operator 
H, (Q 2 ) associated with a hadronic current j. Per- 
turbatively: 



{Hj(Q 2 )}pert = (parton model) 
■ (l + Z™ =1 a n (a s (Q 2 )) n ) 



(8) 



The sum is however formal because of the facto- 
rial growth of the expansion coefficients in large 
orders: 



(a n )uv ~ (-l) n n!6S 



(a n ) 



in 



(9) 
(10) 



where &o is the first coefficient in the /3-function. 
Moreover, the series (0 is associated with inte- 
gration over virtual momenta k 2 3> Q 2 while the 
series (tTU)l is due to k 2 <C Q 2 ■ Assume that the 
perturbative series is asymptotic and introduce 
corresponding power corrections as estimate of 
uncertainties associated with the divergence JHJl 
and (|10|) . It is straightforward to see that the 
corrections are of order A.q CD /Q 2 and Aq CD /Q 4 , 
respectively. 

Now, the central point is that J^J in fact does 
not introduce any uncertainty. The simplest way 
to argue so is to recall that © is Borel summable. 
A more elaborated argumentation exploits 
the fact that at large virtual momenta calcula- 
tions in QCD are reliable because of the asymp- 
totic freedom. 

3.3. Short and long perturbative series 

Let us emphasize again that originally sum 
rules contain short perturbative scries, sec Eq. 
£[J. The 1/Q 2 correction was introduced within 
this context. In fact, the replacement (J7J) is to 
be used in the Born approximation and is not 
suited, say, for the renormalization program. If, 
however, one calculates many terms in the series, 
no 1/Q 2 term can be added since the ultraviolet 
renormalon is calculable. 
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From the point of view of applications, the cru- 
cial question is, how many orders one should cal- 
culate explicitly to avoid adding 1 /Q 2 correction. 
Unfortunately, there exists no theoretical frame- 
work to answer this question and the answer may 
vary from channel to channel. Numerically, the 
gluon condensate is studied by far much better 
than any other observable [201 ■ Perturbatively: 

(0\^(G^) 2 \0) pert = 

const r „ „. , . 
— j- 1 + Sa„< (la)] , 11 
a 4 

where 1/a = Ayy is the ultraviolet cut off (a 
is the lattice spacing) and in the terminology we 
used above, Q 2 — > 1/a 2 . First 10 terms of the 
perturbative expansion were calculated explicitly 
and there is no sign yet that the infrared renor- 
malon sets in. In other words, the 1/Q 2 correc- 
tion is still to be added by hand to describe the 
data. Extrapolation of the coefficients obtained 
to higher n indicates, however, that finally the 
1/Q 2 is absorbed into the perturbative series, as 
is expected theoretically. 

3.4. Comments on literature 

For many years, dominance of soft-field correc- 
tions to the heavy quark-potential at short dis- 
tances was considered obvious, see, e.g. |22 1 - O n 
this background, introduction of a linear term in 
the potential at short distances |6I7I3I24| looked 
heretical. 

Nowadays, it seems to become a common point 
that the soft-field contribution to the heavy-quark 
potential © is negligible and effective linear term 
is there, see, e.g., . Moreover, a careful anal- 
ysis of Ref. [2B1 also supports the guess Q that 
the short distance 1/Q 2 correction is significant 
specifically in one-gluon channel. Once a combi- 
nation of channels is considered which is not con- 
tributed by a one-gluon exchange the soft-field 
contribution might show up 2 . In this sense, the 

1 ln case of bound states of heavy quarks, the effects of re- 
tardation 1221 reduce the correction JFJ by an extra power 
of (Aqcd ■ r) m 

2 After finding this confirmation of 3 6 the authors of Ref. 
1261 summarize their results as being in contradiction with 
these papers. This is a clear misstatement. The author is 



phcnomenological support for the models |(il.'il7| 
has been strengthened recently. 

However, the recent trend |25) is to work with 
long perturbative series which presumably absorb 
the 1/Q 2 corrections. For this reason, direct com- 
parison between and models |6I3| formulated 
within the context of short perturbative series 
is not straightforward. This issue of duality be- 
tween high orders of perturbation theory and the 
quadratic corrections is indeed crucial to compare 
various approaches. 

It is worth mentioning that the recent papers 
|25| concentrate exclusively on the heavy-quark 
potential. No actual 'long enough' perturbative 
series is known explicitly in this case. To substi- 
tute for this lack of knowledge, extrapolations uti- 
lizing infrared renormalons are common. A priori, 
such extrapolations could be correct or wrong . 

The actual check of the models [23, to our 
mind, would be to apply these models in other 
channels, see, e.g., ©. It goes without saying 
that comparison of extrapolations with the ex- 
plicit long series available in case of the gluon 
condensate [20] would be most desirable. 

4. CONFINING FIELDS 

4.1. Monopoles, vortices 

Painstaking search for confining fields in lat- 
tice Yang-Mills theories took many years but was 
crowned with success. Namely the confining fields 
can be viewed either as lattice monopoles, for re- 
view see [2Ej or vortices, for review see More- 
over, the difference between the two languages is 
superficial, see, e.g., J2Hj- Namely, the confining 
field configurations are 2d surfaces (vortices) pop- 
ulated with particles (monopoles). 

Usually monopoles and vortices are treated as 
effective infrared degrees of freedom |28I5| . How- 
ever, they have an ultraviolet facet as well |30I29| . 



thankful to G. Bali and A. Pineda for correspondence on 
this point. 

3 To our knowledge, there is no single example so far of an 
explicit calculation of a perturbative series which would 
exhibit set in of the renormalon-related asymptotic of the 
coefficients. The longest series known 1201 is no exception 
in this sense. For the most recent review of the status of 
renormalons and their possible generalization see 1271 . 
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In particular, for the vortices one finds: 

A « 24(/m)" 2 V4 0, S » 0.5-^ ESI, (12) 

where A is the total area of the vortices, S is non- 
Abelian action associated with the vortices, V4 is 
the total volume. Note that the action diverges 
in the continuum limit of the vanishing lattice 
spacing, a — » 0. 

4.2. Back to the quadratic corrections 

Both monopoles and vortices are extended ge- 
ometrical objects. However, if they are so to 
say projected onto a (quasi)local variable their 
contribution corresponds to the quadratic correc- 
tion considered above [5]. This can be explicitly 
demonstrated in at least two cases. Namely, com- 
bining the two observations cited in (|12fl we get 
for the contribution of the vortices into the gluon 
condensate: 

(0\a s (G^) 2 \0) vortlces ~ — ^j— (Aqcd • a ) 2 ,(13) 

which is nothing else but a quadratic correction 
to the perturbative result. 

Moreover, monopoles and vortices generate lin- 
ear piece in the heavy quark potential at short 
distances: 

lim SV QQ (r) ~ a-r (14) 

where a numerically is very close to the (large- 
distance) string tension. Theoretical reasoning 
for validity of l|14|l was given first in Ref. [21] 
where references to the data can also be found. 

4.3. Duality 

It is worth emphasizing again that the contri- 
bution of the vortices to local observables should 
not be added to perturbation theory but is dual to 
high orders of perturbation theory. In case of the 
gluon condensate this theoretical prediction has 
been already verified. Indeed the difference be- 
tween the full value of the gluon condensate and 
its perturbative value is known explicitly j2HI to 
be of order: 

<0|a s (G^) 2 |0) 4ota , - (0|a s (G^) 2 |0) pert 

" ~ b-QCD ■ (15) 



Note that l|15(l is valid only if more than ten or- 
ders of perturbation theory are subtracted while 
for shorter perturbative series it is the quadratic 
correction that dominates the difference l|15|). 

This duality in evaluating a particular observ- 
able might indicate existence of a dual formula- 
tion of the Yang-Mills theory itself Indeed, 
both vortices and, say, instantons are field con- 
figurations defined in terms of the original Yang- 
Mills fields. Nevertheless instantons are respon- 
sible for terms of order Aq C£) in (|I5|1 and they 
should be added to the perturbative series. While 
the contribution of the vortices is dual to pertur- 
bation theory. The difference can be understood 
if one assumes that the vortices become funda- 
mental variables of a dual formulation of the same 
Yang-Mills theories. 

Indeed, then we should have two alternative 
representations for the same heavy-quark poten- 
tial, in terms of the fields and in terms of the 
strings. In particular the potential at large dis- 
tances is trivial in terms of strings: 

lim Vqq ~ (strings contribution) (16) 

while no close expression is known in terms of 
fields, i.e. perturbatively. Similarly, at short dis- 
tances the potential is trivial in terms of the fields, 
i.e. of the gluon exchange: 

lim Vqq ~ (fields contribution) , (17) 
r— *0 

with no close expression in terms of the strings 
being expected. At intermediate distances, where 
the linear term is visible the potential is calculable 
in two alternative ways, either in terms of strings 
or perturbative series: 

SVqq{t) ~ a ■ r ~ (fields contribution) (18) 
or ~ (strings contribution) . (19) 

Adding the two is not allowed. 

5. DIMENSION-TWO CONDENSATE 

5.1. OPE and condensate of dimension two 

The power corrections in the sum rules QJ start 
with l/Q 4 term, with no 1/Q 2 contribution. The 
reason is that the first gauge-invariant operator, 
that is (G°„) 2 , has dimension d = 4. 
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As a matter of fact, dimension two condensate 
can be introduced \','>'2\ in a gauge invariant fash- 
ion. Indeed, consider first for simplicity magne- 
tostatics and imagine that there exist thin tubes 
carrying magnetic flux. Then the 'vacuum ex- 
pectation value' (0|A 2 |0) cannot vanish. More- 
over, its minimal value is obviously related to the 
magnetic fluxes and gauge invariant. Generaliz- 
ing this logic, one can define gauge invariant con- 
densate of dimension two as the minimal value of 
(0|(A°) 2 |0) in Euclidean space-time. The mean- 
ing of this definition can be scrutinized in great 
detail [3J1 

Moreover, the value (A 2 ) enters the OPE for 
gauge variant quantities, see Ref. [33J and refer- 
ences therein. However, the condensate (A 2 ) m i n 
does not enter the OPE for gauge invariant quan- 
tities. Imagine that one can calculate all terms in 
the OPE for a polarization operator Uj(Q 2 ) in- 
duced by a gauge invariant current j, 

U,{Q 2 ) {partem model){l+E n (A QC D/Q) n ) , (20) 

where Aqcd to a corresponding power is associ- 
ated with matrix elements of gauge invariant op- 
erators and, for simplicity, we omitted perturba- 
tive series. The matrix element ((A^) 2 ) min does 
not appear in the r.h.s. of Eq. 120|) . 

Seemingly, this observation is confused some- 
times for a proof that the condensate of dimension 
two is irrelevant to gauge invariant quantities, say, 
hadronic masses. There is no such proof, however. 
The point is that the expansion (|20")l is asymp- 
totic at best. Thus, it is not a closed expression 
and l|2(J|l cannot determine physical masses in a 
model independent way, as a matter of principle. 
From this point of view, it might be even 'good' 
for (A 2 ) min that it does not enter (|2U[) . Relation 
of the type, m 2 hadron ~ (A 2 ) mm is still perfectly 
possible. 

5.2. (A 2 ) m i n as Higgs condensate 

Thus, it is not surprising at all that there ap- 
pear models which relate (A 2 ) m i n to physical ob- 
servables and these models turn to be successful, 
see, e.g., |%4*) . 

Let me mention specifically an idea that 
{A 2 ) min plays a role of a Higgs field. The prob- 
lem is that the lattice monopoles are defined 



usually j2H] within an Abelian projection of the 
original Yang-Mills fields. In this projection, 
magnetically-charged field plays the role of the 
field describing Cooper pairs in theory of super- 
conductivity. This is not gauge invariant descrip- 
tion however. Moreover, absence of gauge in- 
variant condensate of dimension two might be 
thought of as a general objection to the dual- 
superconductor model of confinement. It was 
conjectured in Ref. [35] that (A 2 ) m i n plays the 
role of Higgs field in the Landau gauge. The 
guess is supported by measurements of variation 
of (A 2 ) m i n inside the confining tube, as function 
of the distance to the central axis of the tube. It 
turns out that the penetration length determined 
in terms of this variation is very close numerically 
to the penetration length in the Maximal Abelian 
gauge determined in terms of the monopole dis- 
tribution. Confirming the idea that (A 2 ) m i n sub- 
stitutes for the Higgs condensate. 

6. CONCLUSIONS 

We have argued that quadratic corrections to 
the sum rules are directly related to the confining 
fields. In terms of the OPE, the quadratic correc- 
tion is a part of the coefficient function in front 
of the unit operator. Thus, to judge on the value 
of the correction one should distinguish models 
with 'short' and 'long' perturbative series kept ex- 
plicit. In other words, the quadratic correction is 
dual to high orders of perturbation theory. As far 
as models are concerned, the model with short- 
distance tachyonic mass (and short perturbative 
series) seems to be most successful phenomeno- 
logically. 
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